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ABSTRACT 

On 9-11 May 1998, BeppoSAX observed the low-luminous Seyfert 1 Galaxy NGC4051 
in a ultra-dim X-ray state. The 2-10 keV flux (1.26 x 10~ 12 erg cm" 2 s -1 ) was about 
20 times fainter than its historical average value, and remained steady along the whole 
observation (~2.3 days). The observed flat spectrum (r ~ 0.8) and intense iron line 
(EWc±600 eV) are best explained assuming that the active nucleus has switched off, 
leaving only a residual reflection component visible. 
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1 INTRODUCTION 

X-ray spectra of Seyfert galaxies above a few keV are dom- 
inated by two components: a power law (hereinafter "pri- 
mary") component emitted very close to the black hole, 
and a secondary (hereinafter "reflection") component aris- 
ing from the reprocessing of the primary radiation by neu- 
tral, optically thick matter (Pounds et al. 1990, Piro et al. 
1990). The reflection component, produced by a Compton 
scattering, forms a broad hump peaking at about 30 keV, 
and an intense iron Ka fluorescent line at 6.4 keV (Light- 
man & White 1988; George & Fabian 1991; Matt, Perola 
& Piro 1991). While in many Seyfert Is at least part of 
the reprocessing occurs in the inner accretion disc, as in- 
dicated by the observed relativistic effects in the iron line 
profile (Tanaka et al. 1995; Nandra et al. 1997), a further re- 
processed component can arise if the nucleus is surrounded 
by any large amount of distant, optically thick and neu- 
tral matter (Ghisellini, Haardt & Matt 1994), such as the 
"torus" expected in the unification scenario (Antonucci & 
Miller 1985; Antonucci 1993) or the dust lanes visible in the 
high-resolution optical images of Seyfert galaxies (Malkan 
et al. 1998; Maiolino et al. 1998). While the reflection com- 
ponent from the accretion disc should lag the primary emis- 
sion by minutes or hours at most, the lag introduced by the 



distant reflector is probably of the order of at least weeks. 
If the active nucleus suddenly switches off, the latter com- 
ponent would then continue echoing the primary compo- 
nent for some time, remaining for a while the only witness 
of the active nucleus past activity. We have observed what 
could be precisely this situation in the nearby (z = 0.0023), 
low-luminosity (< L2-iokcV >5x 10 41 erg s" 1 ) Seyfert 1 
galaxy NGC4051. BeppoSAX observed the source on 9-11 
May 1998 and measured a constant flux along the whole 
observation (~ 2.3 days), which was a factor of five lower 
than the faintest state ever observed (which lasted only 
a few thousands of seconds, Uttley at al. 1998). A con- 
temporaneous and longer EUVE observation suggests that 
such a faint state lasted for at least one week (Fruscione 
et al., in preparation). The spectrum above a few keV is 
completely dominated by a "bare" reflection component 
and implies a luminosity of the primary nuclear compo- 
nent L nuc < 0.14 x 10 41 erg s" 1 . This is the first detection 
of a so dramatic and long-lasting fading in luminosity from 
NGC4051, and provides one of the first direct evidence for 
the presence of large amount of circumnuclear cold, thick 
matter in Seyfert Is environment. 
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NGC4051 was observed by the Italian-Dutch satellite Bep- 
poSAX (Boella et al. 1997a) from 1998 May 9 09:56:44 UTC 
to May 11 18:05:06 UTC. In this paper, data of the imag- 
ing instruments (Low-Energy Concentrator Spectrometer, 
LECS, Parmar et al. 1997; Medium- Energy Concentrator 
Spectrometer, MECS, Boella et al. 1997), and of the col- 
limated Phoswitch Detector System (PDS, Frontera et al. 
1997) are presented. Cleaned and linearized event files have 
been obtained with the reduction software package Saxdas 
(version 1.3.0), using standard screening criteria (see e.g. 
Matt et al. 1997). The PDS data have been further screened 
with a temperature-dependent rise-time threshold, which 
allows a ~ 50% reduction in instrumental background. The 
total net exposure times are 55.9, 69.2 and 33.0 ks for the 
LECS, MECS and PDS, respectively. 

In the following, uncertainties are quoted at 90% level 
of confidence for one interesting parameter (Ax 2 = 2.71); 
energies are in the source rest frame; Ho = 50 km s _1 Mpc -1 
is assumed throughout. 
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Figure 1. MECS and PDS spectra (upper panel) and residuals in 
units of standard deviations (lower panel) when a "bare" Comp- 
ton reflection continuum model + Gaussian line model is applied. 



3 SPECTRAL ANALYSIS 

LECS (MECS) spectra have been extracted in circular re- 
gions of 4 (2) arc minutes and rebinned in order to have at 
least 3 (2) energy channels per resolution element. Back- 
ground spectra from blank sky fields and response ma- 
trices publicly available at the Beppo SAX Science Data 
Center (SDC, September 1997 release) have been employed 
throughout. PDS spectra have been obtained by subtraction 
of the "off-" from the "on-source" intervals. The net count 
rates are (1.12 ±0.05) x 10" 2 s _1 in the LECS (0.1-4 keV), 
(1.09 ± 0.04) x 10" 2 s _1 in the MECS (1.8-10.5 keV) and 
0.123 ± 0.027 s" 1 in the PDS (13-50 keV). No variability 
is observed both in the 0.1-3 and 3-10 keV energy bands, 
contrary to the factor up to 50 typically observed in various 
X-ray energy bands (McHardy et al. 1995; Guainazzi et al. 
1996; Uttley et al. 1998; Fruscione et al., in preparation). 

A fit of the 1.8-10.5 keV MECS spectrum alone with 
a simple power-law model, absorbed through a column 
of neutral Hydrogen with N H = N HGal = 1.3 x 10 20 cm" 2 
(Elvis, Lockman & Wilkes, 1989), yields a very poor fit 
(X 2 = 99/36 de grees of freedom, dof). The main contri- 
bution to the residuals is due to a large emission feature 
at ~ 6.4 keV. Adding a Gaussian profile to the contin- 
uum significantly improves the quality of the fit (Ax 2 = 
46). Best-fit parameters are: photon index P = 0.78±J!;i3, 
line centroid E c = 6.57±o!il keV, line Gaussian dispersion 
a = 0.3 ± 0.3 keV, line equivalent width EW = 1.9±£:§ keV. 
The 2-10 keV flux is ~ 1.4 x 10~ 12 erg cm" 2 s" 2 . Despite the 
approximate nature of the spectral description, three points 
are evident: a) BeppoSAX observed NGC 4051 in an ultra- 
dim state, with a 2-10 keV flux about a factor of 20 fainter 
than its historical average (see Sect. 4); b) the intermedi- 
ate X-ray spectrum is much flatter than in any previously 
reported observation (see e.g.: Fiore et al. 1992; Mihara et 
al. 1994, Guainazzi et al. 1996) and in any other Seyfert 1 
galaxy; c) the iron line EW is several times higher than typi- 
cally observed in Seyfert Is in their "normal" state (Nandra 
& Pounds 1994; Nandra et al. 1997) and in NGC4051 in 
particular (Mihara et al. 1994; Guainazzi et al. 1996). 



These pieces of evidence suggest that the primary nu- 
clear emission has switched off leaving the Compton reflec- 
tion component (along with the iron fluorescent line) as an 
echo of the past activity. We therefore fitted the MECS 
(above 4 keV) and PDS spectra with a "bare" reflection 
component only plus a Gaussian. We used the Xspec model 
pexrav, assuming an inclination of 0°, solar abundances, and 
no high-energy cut-off in the intrinsic power-law; such as- 
sumptions do not substantially affect the following results. 

This model yields a significantly better fit (x 2 = 1.17) 
than the ones where the continuum is fitted with a sim- 
ple power-law (x 2 = 1.38), or with a strongly absorbed 
(N H ~ 3 x 10 24 ) Seyfert-like (r = 1.9) power-law + reflec- 
tion from a plane-parallel infinite slab (x 2 = 1.36), which 
both significantly underestimate the PDS counts below 
30 keV. Table ^ summarizes the best-fit results. The spec- 
tral index of the illuminating (but now invisible) power-law 
(r ~ 1.92, see Figure |l|) turns out to be perfectly consis- 
tent with the typically observed value for NGC4051 (see 
Guainazzi et al. 1996 and references therein). The 2-10 keV 
flux (luminosity) as inferred from the best-fit is (1.26 ± 
0.06) x 10~ 12 erg cm" 2 s" 1 ([2.84 ± 0.14] x 10 40 erg s" 1 ). 
The 90% upper limit on the 1 keV relative normalization 
between the primary and reflected continua is ~ 2%. This 
implies a nuclear flux decrease of at least a factor of 35 in 1.5 
years, if compared with the average measured in the latest 
RXTE long- look pointing (Uttley et al. 1998). Alternatively, 
a column density of at least 1.4 x 10 26 cm -2 is needed, if the 
fading is due to absorption rather than a switch off of the 
primary continuum. 

The line centroid energy is consistent with fluorescence 
from neutral or mildly ionized iron (E c = 6.50±8:io keV), 
with an EW ~ 600 eV. The width of the line is uncon- 
strained, the 90% upper limit on a being an inconclusive 
320 eV. However, some residuals in the blue wing of the 
line may suggest a line blending. Adding a further narrow 
emission line improves the quality of the fit only at a 93.7% 
confidence level. If we assume that an iron Ka/Kp pair is 
produced in the same medium (and therefore with the same 
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Table 1. Best fit parameters ior the May 1998 observation. Errors correspond to A% =2.7. The iron line energy is given in the source 
rest frame. 
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Notes: PO is a power-law, CR is the "bare" Compton reflection continuum, NL is a narrow Gaussian emission line, BL a Gaussian 
emission line where the intrinsic width is left as a free parameter. Photoelectric absorption with Nh = Nn Gal = 1.3 X 10 20 cm -2 has 
been added to all the models above, except for the model # 3, for which Njj = (2.9 ± 0.9) X 10 24 cm -2 . In model # 3 the relative 
normalization between the primary and Compton reflected component has been held fixed to 1. 



intrinsic width), the best fit ratio of the Kp/Ka intensi- 
ties is 0.5±o;4, higher than, albeit not formally inconsistent 
with, the theoretically expected value (~ 11%). Interestingly 
enough, the line flux (I ~ 1.0 x 10 -5 photons cm -2 s -1 ) is 
only a factor of 3 lower than observed by ASCA (Guainazzi 
et al. 1996) or RXTE (Uttley et al. 1998) in NGC 4051 
"normal" states. This shows that a substantial contribution 
to the iron line is likely not to be produced in the imme- 
diate vicinity of the black hole, as previously suggested by 
Guainazzi et al. (1996) and Uttley et al. (1998). 

A prominent excess above the extrapolation of the hard 
X-ray best-fit is evident below 4 keV A detailed spec- 
tral description of the soft X-ray spectrum is beyond the 
scope of this letter, and is deferred to a forthcoming pa- 
per. We only note that if a simple power-law model is 
added to the model CR+BL in Table the soft X-ray 
spectral index turns out to be r so ft = 3.0±o; 2 totally in- 
consistent with the intrinsic spectral index of the primary 
nuclear emission (r 

primary — 1.75±8i|, x 2 = 116.9/93 dof). 
There is no evidence of absorption edges and/or emission 
lines, the 90% upper limit on the optical depths of the O 
vil or O viii photoabsorption edges being 0.34 and 0.13, 
respectively. The 0.5-2 keV (0.1-2 keV) observed flux is 
~ 7.4 (12.1) x 10~ 13 erg cm~ 2 s~\ corresponding to an 
unabsorbed rest frame luminosity of 1.7 ( 7.0) x 10 40 erg s _1 . 



4 THE GLORIOUS PAST OF NGC4051 (AND 
ITS HUMBLE PRESENT) 

NGC 4051 is well-known to exhibit large X-ray flux vari- 
ability on relatively short (~ 10 4 s) timescales (Lawrence et 
al. 1987; Matsuoka et al. 1990; Guainazzi et al. 1996; Uttley 
et al. 1998). BeppoSAX instead observed NGC 4051 in an 
ultra-dim and steady state for the whole ~ 2.3 days of the 
observation. The measured average 2-10 keV flux is more 
than one order of magnitude lower than any historical pub- 
lished flux (with the exception of a poorly constrained EX- 
OS AT 1985 measure, see Figure^). The only comparable ex- 
ample of such a fading is the Narrow Line Emission Galaxy 
NGC 2992, which showed a flux decrease by a factor of ~ 20 
in 16 years (Weaver et al. 1996). This dramatic flux decrease 
is associated with an equally dramatic spectral change. In 
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Figure 2. NGC4051 historical 0.1-2 keV (upper panel) and 2- 
10 keV (lower panel) flux curves, obtained from the high-energy 
X— ray spectra catalog of Malaguti et al. 1994, except for the 1990 
ROSAT/Ginga (Walter et al. 1994), 1991 ROSAT (McHardy et 
al. 1995), ASCA (1993, Mihara et al. 1994; 1994, Guainazzi et al. 
1996) and RXTE long-look (Uttley et al. 1998) observations. The 
error bars represent the 1-er uncertainties on the average fluxes. 

Figure [| we show the 4-50 keV best-fit models as observed 
by ASCA (1994, Guainazzi et al. 1996) and Beppo-SAX. 
In its "normal" state, the NGC4051 2-10 keV spectrum is 
dominated by a continuum with typical intrinsic spectral in- 
dex of F ~ 1.9, to which a Compton reflection component is 
superimposed, flattening the spectrum above ~ 7 — 8 keV. 
This BeppoSAX observation is readily explained if the nu- 
cleus has switched off leaving the reflection component as 
the only witness of the past brightness. This interpretation is 
further supported by a contemporaneous and even longer (8- 
15 May) EUVE (0.02-0.21 keV) monitoring, which caught 
NGC 4051 in a much fainter and constant state then usual 
(Fruscione et al., in preparation). The implied decrease of 
the nuclear luminosity (L nuc < 0.14 x 10 41 erg s _1 ) may be 
explained, if the accretion occurs through a disc, either by a 
proportional (to the luminosity) change in the accretion rate 
or by a smaller change but sufficient to trigger a transition 
of the disc from a radiatively efficient Shakura-Sunayev disc 
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Figure 3. Best-fit model for the 1996 ASCA and 1998 Bep- 
poSAX observations. Solid lines mark the total spectra, while 
dashed lines the spectra components of the ASCA best-fit con- 
tinuum only (for which reflection from a plane-parallel infinite 
slab is assumed). 

to an inefficient advection-dominated flow. Alternatively, it 
may provide support to the idea that thermal-viscous hy- 
drogen ionization disc instabilities play a major role in the 
onset of the AGN phenomenon, which could represent the 
"outburst" phase of accretion from unstable discs (Burderi 
et. al. 1998). 

The present, reflection dominated X-ray spectrum of 
NGC4051 closely resembles those observed in Compton- 
thick Seyfert 2 galaxies (Iwasawa & Comastri 1998; Iwa- 
sawa, Matt & Fabian 1997; Matt et al. 1996, 1997). In those 
cases, the nucleus is likely to be completely obscured rather 
than switched off, a possibility which cannot be formally ex- 
cluded by present data. It requires that a cloud of at least 
Nh > 1.4 x 10 26 cm" 2 has obscured the nuclear region in 
coincidence with the BeppoSAX pointing. It is not easy to 
imagine that such an amount of absorbing matter could be 
dynamically active and remain stable under the combined 
action of the gravitational and radiative pressures in the ac- 
tive nucleus environment. Therefore, the closest similarity 
occurs with our own Galactic Center, where the reflection 
spectrum and the 6.4 keV fluorescent line from the SgrB2 
molecular cloud have been interpreted as the signature of 
past activity in the now quiescent Galactic Center (Sunayev 
et al. 1993; Koyama et al. 1996; Sunayev & Churazov 1998). 

The soft X-rays have undergone a decrease similar to 
the hard X-rays (see Figure ^|). The steepness of the spec- 
trum and the lack of any "warm absorber" imprinting, 
strongly point against a nuclear origin of the soft X-rays. 
We are therefore likely to observe a "remnant" of the faded 
nuclear activity, normally underlying the overwhelming pri- 
mary continuum. The intrinsic 0.1-2 keV luminosity inferred 
by the present data is ~ 7 x 10 40 erg s _1 , a luminosity high 
but not uncommon for a "normal" galaxy. It is therefore 
possible that the dimming of the nucleus has left the galac- 
tic components (supernovae winds, hot halos) as the bulk 
of the observed soft X-rays. Alternatively, the observed soft 
X-rays could have the same origin as the extended (spa- 
tial scale ~ 100 pc) emission observed by the ROSAT HRI 
in NGC 4151 (Morse et al. 1995) and coinciding with the 



thermal emission from a hot (T ~ 10 7 K) and low density 
(n c < 1 cm -3 ) gas, in pressure equilibrium with the Narrow 
Line Region clouds. High energy and spatial resolution ob- 
servations with AXAF or XMM could likely contribute to 
settle this issue. 
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